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Abstract 

To evaluate the weak decays of b-hadrons, the AB — 1 weak effective 
Lagrangian is the foundation. Any new physics beyond the standard 
model (SM) would contribute to the effective Lagrangian through the 
loop integration at the weak scale and evolution from the weak scale 
down to the hadronic scale. In this work we present a systematic analysis 
on the effective Lagrangian which mediates hadronic \AB\ — 1 processes 
in the framework of the minimal flavor violation supersymmetry as well 
as a numerical evaluation of the Wilson coefficients in the effective theory. 
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1 Introduction 

The forthcoming B-factories will make more precise measurements on the rare B-decay processes and those 
measurements would set more strict constraints on the new physics beyond the SM. The purpose to investigate 
B-decays, especially the rare decay modes is to search for traces of new physics and determine its parameter 
space. 

The new physics effects on the rare B processes are intensively discussed in literature. If we believe that the 
SM is only an effective theory and the supersymmetry is more fundamental, measurements of rare B-processes 
will definitely enrich our knowledge in this field. But before we can really pin down any new physics effects, 
we need to carry out a thorough exploration in this area, not only in SM, but also in many plausible models, 
especially the supersymmetric model. 
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The calculation of the rate of inclusive decay B — > X s j is presented by authors of JTJ, |3[ in the two-Higgs 
doublet model (2HDM). The super symmetric effect on B — > X s ^ is discussed in |4|, [B], [(§ and the next-to-leading 
order (NLO) QCD corrections are given in @. The transition b — ► 577 in the supersymmetric extension of 
the standard model is computed in ||. The hadronic B decays^] and CP-violation in those processes[lC] 
have been discussed also. The authors of [11] have discussed possibility of observing supersymmetric effects 
in rare decays B — > X s j and B — > X s e + e~ at the B-factory. Studies on decays B — > (X, K*)l + l~ in SM and 
supersymmetric model have been carried out in |12|. The supersymmetric effects on these processes are very 



interesting and studies on them may shed some light on the general characteristics of the supersymmetric 



model. A relevant review can be found in [ 13 1 . For oscillations of B — Bq (Kq — Kq), calculations have 
been done in the SM and 2HDM. As for the supersymmetric extension of SM, the calculation involving the 
gluino contributions should be re-studied carefully for gluino has a nonzero mass. At the NLO approximation, 
the QCD corrections to the Bq — Bq mixing in the supersymmetry model have been discussed recently. The 
authors of [14, 15] applied the mass-insertion method to estimate QCD corrections to the Bq — Bq mixing. The 



calculations including the gluon-mediated QCD were given in[l6[], and later we have re-derived the formulation 
by including the contribution of gluinos [|l7]] . 

The supersymmetry effects influence the rare B processes in two ways: 

• the Wilson coefficients of the operators existing in the standard model case receive corrections from the 
supersymmetry sector. As far as the four-quark operators are concerned, the supersymmetric contribu- 
tion begins at the order of 0(a s ) along with the SM QCD corrections. 

• when the supersymmetry effects are taken into account, the operator basis is enlarged, new operators 
emerge. 

Generally speaking, theoretical predictions on the inclusive decay rates of B-mesons rest on solid grounds due 
to the fact that these rates can be systematically expanded in powers of QCD [|T^, |l9| , where the leading term 

A 2 

corresponds to the decay width of free b-quark. As the power corrections only start at 0{ Q % ), they affect 

m b 

these rates by at most a few percent. Theoretically, the non-spectator effects of order 16n 2 ( ^ QCD ) could be 
large [20, 21], especially for the charmless decay modes of B-mesons [^, p3| . The main contributions to 



the lifetimes of B-mesons and A& are from the b-quark decays which are thoroughly studied in the framework 
of SM. 

All the theoretical calculations are based on the weak effective Lagrangian which determines the effective 



vertices in the concerned Feynman diagrams[24|. The NLO calculation has been carried out in SM[25, 26], 
but is not complete for the supersymmetric extension. In order to study the supersymmetry effects in those 
low energy processes, one should obtain a complete effective Lagrangian which includes the supersymmetric 
contributions at the order 0(a s ). Here we consider the supersymmetric model with minimal flavor violation 
(MFV), i.e. all flavor transitions occur only in the charged-current sector and are determined by the standard 
Cabibbo-Kobayashi-Maskawa (CKM) mechanism. 

In this work, we perform a complete analysis on the \AB\ = 1 effective Lagrangian, including the current- 
current operators and penguin-induced operators within the framework of the MFV supersymmetry. 

After matching between the full MFV supersymmetric theory and the effective Lagrangian, the Wilson 
coefficients for concerned operators are obtained at the weak scale. Using the recently developed two-loop 



QCD anomalous dimension matrix of flavor changing four-quark operators [27], we discuss the evolution of the 
I AS I = 1 non-leptonic effective Lagrangian from the weak scale down to the hadronic scale. 

The paper is organized as follows. In section 2, we review the minimal flavor violation supersymmetric 
model and give the notations adopted in our analysis. In section 3, the detailed derivations of the effective La- 
grangian at the weak and hadronic scales are made. Then in section 4, we present the numerical results which 
explicitly demonstrate the difference of the Wilson coefficients in MFV supersymmetry and SM. We discuss 
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the numerical results and make a short summary in section 5. Some complicated and tedious expressions are 
collected in the appendices. 



2 The supersymmetry with minimal flavor violation 



Throughout this paper we adopt the notation similar to Ref . ||28f , the expressions of the concerned propagators 
and vertices can be found in the appendix of Ref. p8|] . For convenience, we write down the superpotential 
and relevant mixing matrices. The most general form of the superpotential which does not violate gauge 
invariance and the SM conservation laws is 

W = iieijHjH] + e vj h\ H^R 1 - h d (HlQi - H^Q^D 1 - h T u (H 2 V* J1 Q J 2 - HlQW ■ H) 

Here, the weak SU(2) doublets of quark superfields have been written in the form 




where 1 = 1, 2, 3 are the indices of generations. The Higgs and lepton doublets are denoted by H , H 2 and 
L 1 , respectively. The rest superfields U 1 , D 1 and R 1 are quark superfields of u- and d- types and charged 
leptons in singlets of the weak SU(2). Indices i, j are contracted for the SU(2) group, and hi, h Uj d are the 
Yukawa couplings. 

In order to break the supersymmetry, the soft breaking terms are introduced as 
Csoft = -m 2 m H}*H} - m 2 m H 2 *H 2 - m 2 Ll Li*Lj 

-m 2 Rl R**B* - m 2 Ql Q{*Ql - m^U^U 1 

-mpjD 1 *!) 1 + (miAsAi + m 2 A^A^ 

+m 3 A£A£ + h.c.) + '/>>,,/////; + eijAlhlHmk 1 

-AWAHlQi - H\V lJ Q{)& - Aihi(H 2 V* JI Qi - H 2 Q{)\] 1 + h.c] +■■■ (2) 

where rnr H \, m 2 H 2, m 2 L i, m \i> m Qii m ui an d are t ne parameters of dimension two, while 771,3, m 2, mi 
denote the masses of Xq (a = 1, 2, • • • 8), A^ (i = 1, 2, 3), and Xb, the SU(3) x SU{2) x U(l) gauginos. 
B is a free parameter of dimension one. Af, A^, A d (I = 1, 2, 3) are the soft breaking trilinear couplings of 
scalars. The dots in Eq.|2| stand for flavor-off-diagonal terms (e.g. mQ K Q[*Q{ (5 KI 5 KJ — V KI *V KJ )) that are 
assumed to be negligible. Such terms do occur in our numerical calculation in section 4, but are indeed very 
small. 

Taking into account of the soft breaking terms Eq.(|2|), we can study the phenomenology within the 
supersymmetric extension of the standard model with minimal flavor violation (MFV MSSM). Once the 
favor-off-diagonal soft supersymmetry breaking terms are neglected, the squark mass matrices can be written 
as 2 x 2 matrices for each flavor separately: 

m 2 ^ ( m Qi + m ui + (5 ~ § sin 2 %) cos 2/3m| -m u i (A^ + fj, cot /3) \ 

m u' ~ y -m u i(A I u + fi cot P) m 2 lJl + m 2 7 + | sin 2 % cos 2/3m| J ' 

and 

2 V m Qi + m % + (l + |sin 2 ^ w )cos2/3m| -m d i (A d + /u tan (3) \ 

~ I -rn^j (A d + 11 tan P) m 2 Dl + 777 2 7 - \ sin 2 9 W cos 2/3m| ) ' 
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with m u i, rrifji (I = 1, 2, 3) are the masses of the I-th generation quarks. 

The SM and the MSSM differ in their Higgs sectors. There are four charged scalars, two of them are 
physical massive Higgs bosons and other are massless Goldstones. The mixing matrix can be written as: 

Zh=( Sm( l - C ° S f) (5) 
\ cosp snip J 

with tan/? = ^ and V\,V2 are the vacuum expectation values of the two Higgs scalars. Another matrix 
that we will use is the chargino mixing matrix. The supersymmetric partners of the charged Higgs and 
combine to give two Dirac fermions: xti X2 ■ The two mixing matrices appearing in the Lagrangian are 
defined as 

(Z~) T M C Z + = diag(m Xl , m X2 ), (6) 

where 

M c = 



2ni2 72 92V2 



is the mass matrix of charginos with 52 denoting the gauge coupling of SU(2). In a similar way, Zjj^d 
diagonalize the mass matrices of the up- and down-type squarks respectively: 

Z^m^jZjji = diagimly^rnlji) , 

Z^m^Zjji = diagim^m 2 ^) . (7) 
With those mixing matrices defined above, we can write the interaction vertices as in Ref . pq| . 



3 Matching the coefficients of operators 

As in the SM case, we need to obtain the low-energy effective Lagrangian with five quarks, and while deriving 
it, the heavy supersymmetric degrees of freedom as well as that of SM, including top quark, W-bosons, charged 
Higgs bosons and the supersymmetric partners of the standard particles are integrated out. In this work we 
only retain the operators up to dimension six. In this approximation the effective Lagrangian for \AB\ = 1 
reads 

C eff = -^-V: s V tb [Y: (ct{n)Qi + CT(/i)4 c ) + E c ^)Qf\ (8) 



where Gf is the Fermi coupling constant, Cf(fi), Cf(fi), Cj(fi) (i = 1, 2, • • • , 10; j = 1, • • • , 5) are the 
Wilson coefficients evaluated at the scale fi; Vtb and Vt s are the matrix elements of the CKM matrix. Making 
the effective Lagrangian close under the QCD renormalization, we include the penguin operator beside 
those four-quark operators. 

As commonly adopted in literature, we classify the operators as the current-current operators which are 
originally induced by the tree level W-exchange interaction and one-loop 'box' diagrams, and the "penguin- 
" induced operators. Both of them would undergo QCD corrections and receive contributions from supersym- 
metric particles via loops. 
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The current-current operators are written as|27] 



Ql 

Ql 
Q c 3 
Ql 
Ql 
Ql 
Q c 7 
Qs 
Ql 

Qto = ( 



oj-bi 



) (c/37^ 



S Q£ UJ — Cqi 



) {cpa^uj-bo^ , 

^ (cp^u+bp) , 
S a U^Cp^(cpU + b a J , 
(cpLO+bpj , 



S qi UJ — Cfj 



(9) 



where = §[7^,7^ and lo± = 1 J 5 . In the standard model, there are only two such operators i.e. Q\, Q2, 
18 new operators are induced when super symmetry takes part in the game. There are other ten current- 
current operators which are simply obtained by interchanging uj± <-> uj t in Qf, i.e. Qf = Q1(uj± <-» Due 
to the small CKM entries for the u-analog operators in the |AJ3| = 1 effective Lagrangian, V* s V u b « V* s V c b 
the u-quark analogs of the effective Lagrangian in Eq.|| can be neglected. The basis of the penguin-induced 
operators consists of [29, 3C, 31, 52, 33] 



Ql = (s Q 7/^-&a) Y 

g 

Ql = {sal^-bpj Yl (<7/37 M ^-<?a) , 
Ql = (saln^-ba) Y (9/37 M ^+^) » 



Ql = (aaJnU-bp) Y ( 



QpT^+Qa 



Ql 



(4tt) 2 



(10) 



.„-sg s G ■ <j{m s u- + m h u + )b , 

r 

with q = u, d, c, s, b, G^ v = G a ^ v T a denotes the gluon field strength tensor, G°^ v = d^-d^+gsf^G^Gl, 
andC'-rr G^a'"'. 

In the following sections, we will derive the Wilson coefficients for the current-current operators in Eq.^ 
and penguin-induced operators in EqJlO|. 
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3.1 The difference of quark field normalization in the full and effective theories while 
taking the MS scheme 

In order to systematically investigate all corrections of the supersymmetry QCD to the vertex duW including 
the self-energies of external quark legs, Ciuchini et.al used the on-shell scheme to subtract the divergence of 
the quark fields and the duW vertex that originates from supersymmetric partners, whereas the divergence 
originating from the SM sector is still subtracted out in the MS-scheme@. In this work, we will employ the 
MS scheme throughout and show that the results are qualitatively consistent with theirs. 

It is noted that there is a difference of the normalization of the quark fields in the full supersymmetric 
theory and the effective theory while taking the MS scheme. As a matter of fact, for the self-energy which 
determines the renormalization of the wave functions, the Feynman diagrams in the effective theory are the 
same as the standard model part of the full theory, therefore there is not normalization ambiguity in the SM 
case. When the supersymmetry sector is included, there exists an extra term from the supersymmetry sector 
and the normalizations of the external quark fields are not the same after renormalization in the MS scheme, 
as 

* 1 



where AZ^' = AZ e ?f and the superscript "eff" denotes the quantities in the effective theory. Thus one 
cannot simply match the vertex-induced Lagrangian in the full and effective theories because the external 
quark fields have different normalizations. That is understood that in both the full and effective theories, 
the SM quarks exist, but the supersymmetry particles (squarks and gluinos) only exist in the full theory, 
but are integrated out to produce the effective Lagrangian. The difference AZf ul1 - AZ e ^ = AZ{f^ SM 
is a finite renormalization effect which should be included when we match the Lagrangian in the full and 
effective theories. Namely, when we match the Lagrangians, we not only consider the contributions from 
the vertices, but also include this normalization difference. We find that the large logarithms which exist 
in the vertex contributions and this normalization difference would cancel each other exactly and then the 
decoupling theorem is obvious. 

The difference manifests as a finite renormalization contribution to the AB = 1 effective Lagrangian, 
which will be expressed explicitly in the following computation. 

The supersymmetric contributions to the self-energy are 



( , o x.j x.Anx. r xz hxx. — x.x, , \nx„ T N 

(p) = -i^-Cpl A + lnx M + - + *— + — s *- >-2 " ( Z bi Z vtfi 

4vr [V 2 x--x fjI x--x fiI (x- - x &1 ) 2 J\ u u 



OLJ- 



+Z^Z^ + )+2m ^ ^ lnX&/ (z$ZlW + Zl i I Z%?u J+ )\ , 
( / o X-r X- r lnx- r x 2 . lnx- — x-x . , Inx . , 

a. „ / , , O D 1 . D 1 D 1 9 9 9 D 1 nl 



iH d (p)=-i-lC F {[A + \nx f , + - + - ^— + 



47r V 2 x- — x. r x-—x. r (x~—x~ T ) 

( x 9 nl 9 nl v 9 £>I ' 




+44>+) + 2m ' 9 X6{il l X [ x lnX6{ \ z 5* z >- + z ^ z > 

9 b 1 

t 

where x,, = ^W- , x; = ^j- . In the expressions the first term is the correction to the wave function of the 

'"w W 

quarks whereas the second term corresponds to the supersymmetric contributions to the quark masses. In 



G 



Eq.ll, if we complete the renormalization for the quark fields in the MS-scheme, the divergent part in the 



renormalization multiplier which should be subtracted is 

8Z±(MS) = —C F A . (12) 
H Air 

When matching the full and effective theories in MS scheme, there is a difference of the renormalized quark 
fields, so that we need to derive the extra contribution. 

In literature [pj], the on-mass-shell renormalization [44| is adopted, by which after renormalization the 



residue of the self-energy should be one at the physical mass, so that the normalization ambiguity does not 



exist. However, in the MS scheme, there is no such a requirement. From Eq.ll, we have the difference of the 
normalization of the quark fields in the full supersymmetric theory and the effective theory while taking the 
MS scheme as 

/ o x.j x- r lnx_ r x 2 Inx. — x~x. . Inx, , 

Ols „ /, O „1„- „1„-^/ U 1 U 1 U 1 9 3 3 V l uf 



4f \ 2 u u \x.—x, T x~—x. r (x-—x. r ) 

I 9 u! 9 Jjl y 9 If 1 ' 
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, o x.j x„ T lnx„, x 2 lnx. — x.x . , Inx . , N 

AZ dI = £Wlnx M + | + £ Z%Z% (-^- + -*-^L - 9 ' '% * )) , (13) 

a 4"7T \ ^2^^^ \X-—X- T X-—X-, [X.—X-jY 'I 

x % 9 D 1 9 D 1 9 

i i % 

When we match the Lagrangians in the full and effective theories, we not only consider the vertex-induced 
contributions, but also need to include this finite renormalization contribution. Thus taking this normalization 
difference into account, we obtain the Wilson coefficients of the current-current operators in the full and 
effective theories. 

3.2 The Wilson coefficients of current-current operators 

The Wilson coefficients Cf, Cf (i = 1, 2, • • • , 10) in Eq]8| can be determined by the requirement that the 
amplitude Af u n in the full theory is equal to the corresponding amplitude in the effective theory at the weak 
scale 



A full = A eff = ^-V^V cb f^{cfQt + CtQt - [AZ- + l -AZ b ~ + \aZj]qi). (14) 



The QCD induced one-loop Feynman diagrams responsible for AB = 1 effective Lagrangian in the minimal 
flavor violation supersymmetric theory and effective theory are drawn in Fig.|l| and Fig.Q respectively. The 
last term of Eq]l4] originates from the difference of quark field normalization in the full and effective theories. 
Considering the QCD corrections of current-current operators, we can extract the Wilson coefficients Cf , Cf 
at the //w scale. The expressions for those coefficients are presented in appendix [A|, and one can notice that 
all resultant Cf and Cf are free of infrared divergence. The first terms of Cf (fj-w) and C^i^w) are the SM 
contribution whereas the other terms are due to the supersymmetry contributions. Other non-zero Wilson 
coefficients Cf (nw), Cj (/J-w) (i = 3,4, ...,10, j = 1,2,..., 10) all originate from contributions of scalar quarks 
and gluino. Provided rrig = = m SUSY , we have 

a 11 

C^(n w ) ~ C^ SM (nw) = 1 - Tl( ln ^ + 



4vr p 6 
3a s 11, 



Cf{nw) ~ C% (jiw) = -r^( ln ^ + "7-)> as m SUSY > fj, w , (15) 
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which recovers the SM result. When rrig — ► oo, but m^y = rrifji (I = 1, 2, 3; i = 1, 2) and = 
mjyi {I = 1, 2, 3; i = 1, 2) remain finite, we also have Ci 2 {^w) ~ Ci 2sM (fJ-w)- From the equations given 
in the appendix, one observes that in C^i^w) if all the mass parameters of the supersymmetric particles 
tend to infinity or they are the same, i.e. the supersymmetric particles are degenerate in mass (it is exactly 
the mSUGRA case) all the x-value related terms cancel each other. Here we take into account the fact 
Zjj lt2 ~ -2^1,2,3 ~ I in the mSUGRA scenario. In other coefficients Cf, Cf etc. the logarithm-related terms 
are automatically suppressed to zero when the supersymmetric particles are very heavy (see appendix A for 
details) . Those results indicate the decoupling of the supersymmetric sector as the supersymmetric partners 
turn to be very heavy. 

In order to give a complete |AS| = 1 non-leptonic effective Lagrangian with five quarks, we should include 
the contributions of penguin diagrams. In our present work, we only consider the gluon-penguin. 



3.3 The Wilson coefficients of penguin-induced operators 

In this section, we derive the Wilson coefficients for the penguin-induced dimension six operators. The basis 
for the penguin-induced operators is given in EgjlOj. The one-loop Feynman diagrams for the penguin-induced 
operators are drawn in Fig.||. The obtained coefficients Cf(/iw) read 



CfOw) 
CfOw) 
Cf(/xw) 



4ir 
4ir 

4ir 

a s 
4ir 



- lnx u 
9 



I 1 

-h(x t , x H , , x Xj ) + - 



lnx^ + E(xt, x H , x^ 3 , x Xj ) 



- lnx u 
9 



I 1 

-h(x t , x H ,x f)3 , x Xj ) + - 



lnx^ + E{xt, x H , x^ 3 , x Xj ) 



where 



E(xt, x H , x^ 3 , x Xj ) 



r 18x 4 - llxf - 
12(1 -x t ) 



+ 



(-4 + lGx t - 9x 2 )lnxt 



6(1 -xj) 4 



1 ^(2xtX 3 H — 3x 2 x 2 H )(\nxt — lnx } 



tan 2 



6(x, 



+ 



16x t x 2 — 29x?x„ + 7x 



x t 



x?..(lnx- 3 — In x 



Xj- 



lis 



Xj 



36(x H — 



+ 



18(x 



Xj X (j3 



(16) 



(17) 



The Wilson coefficient of Q§ is 



3xt Inxt 2 + 5xt 



+ 



4(1 -x t ) 



1 



(1 " xt) 



+ 



x t x 2 H (lnx t - lnx H ) 3x H x t 



2{ Xl 



x t ) 



+ 



xt 



4(x H - x t y 



tan 2 j3 



x 2 x 2 H (Inxt — lnx H ) 2x t x 2 H + 5x 2 x H 
4(x H - x t ) 

x?,.X-, (lnx_ 
Xj u? v u 



lnx 



Xj, 



24(x H - x t ) 3 
2x 2 _. + 5x Xj x,. 



Xj 



X 



2(x 



Xj 



x- 3 ) 4 

TTO / 



+ 



12 x 



Xj 



X Q3 j 
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E 



... r^Xj^t/S ( mx ^3 



lnx 



(X 



+ 



X ij3 



Xj 



2 x 



Xj 



(18) 



The first terms of Eq.17 and Eq.18 are the SM contributions! 34, 



% 37 



the second and the third 

terms are the contributions from charged Higgs and supersymmetric particles respectively 0]. In the above 
expression, we only keep the contribution of the up-type scalar quarks of the third generation, and that from 
other squarks are ignored for their large masses. The matrix .A, B are written as 

_ mt 



A ij 



-Z\tZ- 



+ 



\/2mw sin 



B 13 



V^raw cos j3 



mt 



\/2m\v sin f5 



?2i 



(19) 



with Zfj 3 , Z + are the mixing matrices of scalar top quarks and charginos respectively. Because we are 
working in the framework of supersymmetric extension of SM, the effective vertex b — > sg must take in the 
contribution of supersymmetric particles, thus the penguin-induced Lagrangian is somewhat different from 
that within SM. It is noticed that the mn related terms i.e. the last two terms in Eq.17 and Eq.18 tend to 
zero when mu — > oo and as well as m SUSY — ► oo. 

When we only keep the Yukawa coupling of the top-quark to the Higgs boson and assuming that the 
masses of the scalar quarks except top scalar quarks are highly degenerate, and the weak eigenstates are 
the eigenstates of the masses, i.e. Zjj = Zjj = 1, we find that only the coefficients of Q\, Q 2 and the five 
penguin-induced operators are not zero, whereas the Wilson coefficients of other operators vanish. This is 
exactly the operator basis existing in the standard model. In the MS-scheme, the Wilson coefficients are 
simplified as 

9 i X 2 

a, ^ r xtlnx^ 



w 



47T V 



11 

6~ 



In a; 



+ 



{x -x-){x -x-f (x a -x,)(x 



X 



x 2 . In x . 



x. 



'{x 



+ 



x 2 In x - 

Q Q 



x-(21nx- + 1) 



i x s - x /i)( x . 



C 2 c (/i W )=(l-| £ (ln^ + H ) ) 



2vr 



x - 
Q 



x - In x _ 

Q Q 



{x-x 

+ 



+ 



X ){x 



X. 



x 



i In x ■ 



+ x In x- 

Q 9 



X, 



(x s 



X-) 2 



(20) 



the Wilson coefficients of the penguin-induced operators remain unchanged and the other Wilson coefficients 
vanish. 



3.4 The evolution of the Wilson coefficients 

Many flavor-changing processes occur at the hadronic scale with Hhadron *C fJ-w- Generally, first a few 
terms in the perturbative expansion of the amplitudes are sufficient when the renormalization scale \i is 
close to /J-hadron rather than to jU-yy|p^|. The Wilson coefficients Cf, Cf, Cf at Hhadron are obtained from 
C-(/iw), Cf(^w)) Cf(fJ'w) with help of the Renormalization Group Equations (RGEs) evolution. If we define 
a 1 x 24 matrix as 

C~* ( r~ic s~tc f*ic /ntc f~ic f~ic s~iP s~iP s~iy riV\ /oi \ 

— I L/ 1? 2 , • • • , O 10 , 1 , L/ 2 , • • • , L/ 10 , L/ x , U 2 , Lyg , U 4 I , 

and the RGEs for the Wilson coefficients are 

MT"C(/i) = C( M )7(/x) . (22) 
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Here, 7 is the anomalous dimension matrix which has the following form in the perturbative expansion 



= ^7™ + (^j (23) 

Adopting the Naive Dimensional Regularization— MS (NDR-MS) scheme, Buras et.al have given the anomalous 
dimension matrix up to two-loop order in the basis of Eq.|| and Eq|l(]|27], 39]. 



With the RGEs (|22| ) and the Wilson coefficients at the weak scale as the initial condition, we can derive 
the effective Lagrangian of five quarks at the hadronic scale. 

4 Numerical results 

Indeed, there are too many free parameters in the minimal supersymmetric extension of SM (MSSM). In order 
to reduce the number of free parameters, we assume that the MSSM is a low-energy effective theory of a more 
fundamental theory which exists at a higher scale, such as the grand unification scale or the Planck scale. A 



realization of this idea is the minimal supergravity (mSUGRA), which is fully specified by five parameters [40] 

mo, mi, Aq, tan/3, sgn(fi). 

2 

Here mo, mi and Aq are the universal scalar quark mass, gaugino mass and trilinear scalar coupling. They are 

2 

assumed to arise through supersymmetry breaking in a hidden-sector at the GUT scale hgut — 2 x 10 16 GeV. In 
our numerical calculation, to maintain consistency of the theory and the up-to-date experimental observation, 
when we obtain the numerical value of the Higgs mass in the mSUGRA model with the five parameters, we 
include all one-loop effects in the Higgs potential [pi]]. Moreover we also employ the two-loop RGEs[42] with 



one- loop threshold corrections [O, 43] as the energy scale runs down from the mSUGRA scale to the lower 



weak scale. In the framework of minimal supergravity, the unification assumptions are expressed as 

A{ = A I d = Al = A , (24) 
B = A -1, 

22222222 
m H i = m H 2 = m L i = m Rl = m Ql = m ul = m Dl = m , 

777-1 = m 2 = m 3 = m l ■ 

2 

For the SM parameters, we have mj, = 5GeV, mt = 174GeV, m\v = 80.23GeV, a e (m\y) = j^g> a,(mw) = 
0.12 at the weak scale. Taking above values, we find the SM prediction for Wilson Coefficients as 
CT gM (m&) = —0.295, C2 SM (m,b) = 1.110, Ci SM ( m b) = 0.014. In our numerical calculations of the super- 
symmetry corrections to those Wilson coefficients C?(m&) (i = 1, 10), Cf(nib) (i = 1, 10), and Cf(rrib), 
we always set Aq = and sgn{n) = +. Even though other Wilson coefficients also get nonzero contributions 
from the supersymmetric sector, our discussions mainly focus at the dependence of C^m^), C^m^), C\{m\,) 
on the supersymmetric parameters because they play more significant roles in low energy phenomenology. In 
Fig.|I] (a), (b) and (c), we plot the ratios between supersymmetry corrections to Cf (m&), C^m^), Cf (m&) and 

their SM prediction values versus the parameter mi with mo = 200GeV and tan 3 = 2 or 20. The dependence 

2 

of those ratios on mo (mi = 300GeV and tan/3 = 2 or 20) is plotted in Fig.||. 

5 Discussions and Conclusion 

The |AJB| = 1 non-leptonic effective Lagrangian has been considered in the minimal flavor violating super- 
symmetry scenario. The supersymmetry contributions affect the effective Lagrangian via two aspects: 
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• new current-current operators emerge beside those 'old' operators in the SM case; 

• for the 'old' operators, the Wilson coefficients at the weak scale are modified by the supersymmetric 
contributions. 

Now let us briefly discuss our observation of the numerical results. 

In Fig.|| (a) and (b), the two lines (solid and dash) corresponding to tan/3 = 2 and tan/3 = 20 differ 
from each other more obviously. The supersymmetry corrections to the Wilson coefficients Cf 2 ( 771 b) are 
relatively large when the supersymmetry particles have masses of the same order of electroweak energy scale, 
for example, mi = 300GeV, uiq = 200GeV, the supersymmetry corrections to Cf 2 ( m b) can reach about 8%. 

2 ' 

When the masses of the supersymmetric particles become very large, the supersymmetry corrections turn to 
zero due to the decoupling theorem. In Fig.|I](c), the two lines corresponding to tan/3 = 2 and tan/3 = 20 
almost overlap on each other. It is noted that at the left part of Fig. 5 (c) as mi < 2 TeV, there is a sharp 

peak at the dependence of C?(m&) on mi. It drops very fast as mi is away from this region, the resonance is 

2 2 

due to an almost degeneracy of the mass of chargino and the mass of stop and it leads to an obvious deviation 
of the Cf(mb) value from the prediction of SM (0.014). As mi further increases, Cf tends to the predictive 
value of SM. This mass resonance only occurs for the coefficients of the penguin-induced operators, but not 
for the current-current-quark operators, so that the dependence of Cf(m&), Coimb) on mi does not change 

2 

drastically. All these are consistent with the common sense which is familiar to us even before the calculations 
are done. 

In FigJE], we plot the dependence of Cf (m^), C^ivdb) and Cf on mo, with Aq = 0, sgn(n) = +1, mi = 
200 GeV. Similar to the case of Fig]||, the dependence of supersymmetry corrections to Cf 2 (^fe) on mo is 
remarkable and the values obviously deviate from the SM prediction when mo take smaller values. When 
mo ~ 200GeV, Cf 2 { m b) deviate from the SM prediction by 8%, and when mo > 5 TeV, the deviation tends 
to zero. The dependence of Cf(mb) on mo is in analog to its dependence on mi. When mo < 2 TeV, there 
is an obvious peak which damps steeply, the reason is still due to the mass degeneracy of heavy chargino 
and stop. When mo turns very large, Cf (m&) approaches the SM prediction. By contrast, the dependence of 
Cf( m b)) C^imb) on mo is more smooth, because there is no mass resonance effect in this situation. 

In above discussions, we only list the Wilson coefficients of a few typical operators, the Wilson coefficients 
of other operators are in analogy. For simplicity, in all the numerical calculations, we always adopt Aq = 
0, sgn(fi) = +1, but this convention is not necessary. If we dismiss this assumption and let Aq / 0, sgn(fi) = 
±1, the parameter space is enlarged, but the qualitative conclusion remains the same. 

Here we briefly discuss the decoupling of the supersymmetry particles as the concerned energy scale turns 
to infinity. When we match the effective theory to the full theory which is the standard model at the weak 
scale, the quark and gluon fields exist in both theories, thus the renormalized quark fields are the same in 
both theories. However, as the full theory is the supersymmetric extension of SM, scalar quarks and gluinos 
exist only in the full theory, but do not in the effective one. Thus when we match the two theories at a certain 
scale, disappearance of such supersymmetric partners in the effective theory would result in a difference of 
the normalization of the quark fields in the full and effective theories. Actually, this is a finite renormalization 
contribution of the self-energy. Ignoring this effect would lead to a situation that heavy supersymmetry 
particles do not decouple. This normalization difference would also affect evaluation of the vertex-induced 
contribution. When we match the full and effective theories we need to take this normalization difference into 
account, then the large logarithms emerging from the self-energy exactly cancel out that from the vertices, 
thus the decoupling is obvious. On other side, only in the MS scheme, there exists the difference of the 
normalization of the renormalized quark fields, but it does not exist in the on-mass shell renormalization 
scheme. The reason is that requiring the external quarks to be on their physical mass shell and the residue 
for the self-energy to be one can serve as an additional condition by which the normalization of the quark 
fields are the same in the full and effective theories. A direct consequence of the correct renormalization 
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scheme is the decoupling of the supersymmetry sector as the supersymmetry particles are too heavy. If we 
include the normalization difference as we take the MS scheme, the result coincides with that in the on-shell 
renormalization. 

Because we have carefully considered the normalization difference, in our final expressions, one can imme- 
diately observe that the supersymmetry particles would decouple if their mass scale turns to infinity That is 
consistent with the common sense. 

From our discussion and numerical results, one can note that in general, the correction of supersymmetry 
to the Wilson coefficients can be as large as 8% when the supersymmetry particles are not very heavy. It is 
also noted that the |A5| = 1 effective Lagrangian at the weak scale induces an electromagnetic dipole and 
a chromo-dipole. Through the evolution down to the hadronic scale, the super symmetric contribution would 
result in substantial values for the dipoles at low energies which may not be negligible for phenomenology 
when g and Q are relative light. 

With collection of large amount of data at the B-factory and elsewhere main laboratories in the world, 
the measurements on the rare processes b — > 57, b — > suu etc. would set more rigorous constraint on the 
parameter space of the super symmetric model or we can, as expected, find some evidence of existence of 
supersymmetric particles. 
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A The Wilson coefficients of current-current operators at the weak scale 



After matching at the weak scale, the Wilson coefficients for the current-current operators are written as 
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Here, Z ± are the mixing matrices of charginos. 
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Figure 3: The one-loop diagrams for calculating the penguin-induced four-quark operators 



19 



10 




-2 1 1 ' ' ' ' 1 1 ' ' 1 

500 1000 1500 2000 2500 3000 3500 4000 4500 5000 

™ 1 (GeV) 

2 

10 | 1 1 1 1 1 1 1 1 1 1 




- 



-2 1 ' ' 1 1 ' ' ' 1 1 1 

500 1000 1500 2000 2500 3000 3500 4000 4500 5000 

m i (GeV) 

100 ! 1 1 1 1 1 1 1 1 1 1 




500 1000 1500 2000 2500 3000 3500 4000 4500 5000 

m i (GeV) 

Figure 4: The relative super symmetry corrections (the supersymmetric corrections/ the SM predictions) to 
the Wilson coefficients at the scale versus mi with tan/3 = 2 (Solid-lines), tan/3 = 20 (Dash-Lines). The 
other parameters are set as tuq = 200GeV, A$ = 0, sgn(ji) = +. 
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Figure 5: The relative supersymmetry corrections to the Wilson Coefficients at the m& scale versus mo with 
tan/3 = 2 (Solid-lines), tan/3 = 20 (Dash-Lines). The other parameters are set as mi = 300GeV, Aq = 
0, sgn(fi) = +. 
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